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Fermionic functional renormalization group for 
first-order phase transitions: a mean-field model 

Abstract. First-order phase transitions in many-fermion systems are not detected in 
the susceptibihty analysis of common rcnormalization-group (RG) approaches. Here we 
introduce a counterterm technique within the functional rcnormalization-group (fRG) 
formalism which allows access to all stable and metastable configurations. It becomes 
possible to study symmetry-broken states which occur through first-order transitions 
as well as hysteresis phenomena. For continuous transitions, the standard results are 
reproduced. As an example, we study discrete-symmetry breaking in a mean-field 
model for a commensurate charge-density wave. An additional benefit of the approach 
is that away from the critical temperature for the breaking of discrete symmetries large 
interactions can be avoided at all RG scales. 
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1. Introduction 

Renormalization group (RG) methods are powerful systematic tools to investigate 
the emergence of the low-energy physics of many-particle systems. In the Wilsonian 
RG approach [H E], [3] for a many-fermion system, the quasiparticle excitations are 
integrated out successively with an energy scale decreasing toward the Fermi energy. 
This generates an energy-scale dependent effective interaction from which the leading 
interaction processes and correlations can be read off. Although most implementations 
involve a number of approximations, the flow of the vertex functions of the effective 
theory can be derived from an exact RG equation for a generating functional [U [5]. 
This powerful equation sparks the idea that a much more controlled approach to the 
low-energy phase than obtained previously is indeed possible. 

In many applications for correlated electron systems one is primarily interested 
in changes of the ground state toward symmetry-breaking states. If there is such a 
symmetry change of the ground state, a certain bilinear Bij in the fermion operators 
c^^ \ c^p with some quantum numbers i and j takes a nonzero anomalous (i.e. violating 
the orginal symmetry of the action) expectation value at a low RG scale. For a second 
order transition, this occurs continuously as RG scale or temperature are lowered. Then 
the tendency toward symmetry breaking is signaled in the RG flow by a divergence 
of a certain component of the effective interaction at a nonzero energy scale Ac. In 
correspondence with this, the zero-frequency susceptibility for the coupling to the 
bilinear operator Bij diverges in power-law form at the same critical scale Ac. If we 
view the expectation value {B^j) as a bosonic field cf) (e.g. by a Hubbard-Stratonovich 
transformation), the inverse susceptibility appears in the quadratic part of the action 
for the 0-field. The symmetry breaking occurs at the temperature or energy scale where 
this quadratic part changes its sign, i.e. where the effective potential for reverses its 
curvature at = from positive to negative, such that the potential minimum occurs 
at a nonzero value of (f). When the energy scale or temperature is lowered further, 
this minimum moves further away from zero and finally saturates unless collective 
fluctuations move it back to zero and restore the symmetric state (see, e.g., [6]). 

As indicated in this description, the passage from the symmetric state at high 
energy scales or temperatures to the symmetry-broken state at low energies is most 
easily described in a bosonic language, and such methods are being developed ([HinilH])- 
However, in many interesting examples in many-fermion systems it is not clear from the 
outset which type of boson or which fermion bilinear acquires a nonzero expectation 
value. Hence, a purely fermionic description is a less biased approach, as it allows 
the system to develop its intrinsic correlations without any prejudice of grouping 
certain fermions in pairs and treating these in a saddle-point approach. Recently, 
some of us have presented a purely fermionic functional RG method based on [9] 
which allows one to continuously flow into the symmetry-broken state, requiring only 
an infinitesimal symmetry-breaking in the initial conditions. This approach has been 
shown to work correctly for models where mean- field solutions are exact, like the reduced 
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Bardeen-Cooper-Schrieffer (BCS) model [10] or a similar charge-density-wave model 
with discrete-symmetry breaking [111 112] • Extensions for more general models will 
follow. 

Yet, these new fRG extensions face two problems. First, there is still an almost- 
divergence of the interactions at the critical scale. Large values of the effective 
interaction are problematic for the extension to more general models, as they might 
cause a breakdown of the perturbative approach. A true divergence is avoided by the 
small explicit symmetry breaking put into the initial conditions. But in order to stay 
close to the symmetric model which one typically wants to study, some components of 
the interaction in the symmetry-breaking channel still reach large values whose feedback 
on other channels needs to be understood. While this problem may be solved either by a 
compromise with an adequate size of the initial symmetry breaking or by understanding 
the potentially rich and interesting physics of the effect of these large interactions on 
the flow, a second problem is more fundamental. The currently available fRG methods 
only test the stability of the symmetric phase locally around the origin in the space of 
anomalous expectation values. If there is no divergent susceptibility or no curvature 
change in the effective potential of a bosonic field, the fRG flow will remain finite 
and the anomalous expectation values will vanish continuously when the small initial 
symmetry-breaking is sent to zero. This in turn would be interpreted as stability of the 
symmetric phase and absence of spontaneous symmetry breaking. On the other hand, 
simple mean-field treatments of many models indicate first-order transitions into the 
symmetry-broken state. This occurs for example in the two-dimensional t-t' Hubbard 
model near half filling for the (commensurate) antiferromagnetic spin-density wave 
transition [131 [13] and in a model with a particular forward scattering for the transition 
into a state with broken fourfold symmetry of the Fermi surface [15]. The fRG flow for 
these parameters in the almost symmetric state finds an enhanced but not a diverging 
susceptibility, hence missing the symmetry change of the ground state. The fRG method 
simply does not see the thermodynamic-potential minimum for nonzero order parameter 
if it is separated from the symmetric minimum by a finite energy barrier, as in the case 
of certain first-order transitions. What is needed is a method that, intuitively described, 
starts the flow further out in the symmetry-broken range and which then flows into the 
correct minimum corresponding to a spontaneously broken state of a model with a fully 
symmetric action. 

In this work, we present a counterterm extension of the fRG method which is 
able to find symmetry-broken states not captured by the treatment of [HI |12] . These 
would occur via first-order transitions upon varying temperature, chemical potential, or 
infrared cutoff. The procedure invokes a flow differing slightly from Wilson's momentum- 
shell RG. The propagator of the fermionic excitations is increased continuously from 
zero to its full strength with the same rate for all excitation energies. This flow, also 
dubbed interaction flow, has been tested previously for the 2D Hubbard model without 
symmetry breaking and gives analogous results to those obtained with momentum-shell 
or temperature-flow schemes [16]. It can be regarded as an infinitely smooth cutoff, i.e. 
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the opposite extreme to a sharp cutoff, for which all modes above a running energy scale 
A are included with full weight and all modes below A do not contribute. 

Counterterms have been employed in the context of correlated fermions before, 
notably in calculations for the symmetric state of interacting electron systems. The 
early work [T7j by Feldman and Trubowitz considers systems with spherical free Fermi 
surfaces. Unphysical infrared divergencies of perturbation theory Feynman graphs are 
removed by adjusting the chemical potential on the propagator lines to the interacting 
Fermi surface, countering this by an appropriate addition to the self-energy. The 
addition is determined by renormalization-group methods. This procedure can be 
viewed as a shift from the chemical potential into the interaction. Feldman, Salmhofer, 
and Trubowitz [18] employ it for more general Fermi surfaces which necessitate a 
counterterm that is a function in momentum space. They identify and find bounds 
for the most singular contributions to this couterterm function under slightly more 
restrictive conditions [19]. The authors prove regularity properties of the counterterm 
function in [20] , which they use in [21] to study the invertibility of the mapping induced 
by the counterterm which transforms the interacting to the corresponding free Fermi 
surface. The analysis of [21] is restricted to space-inversion-symmetric dispersions. 

Counterterms have also been used in the study of symmetry breaking. Feldman and 
Trubowitz [22] is devoted to flows with symmetry-breaking in the Cooper channel for a 
spherically symmetric electron-phonon system. Counterterms are again used to remove 
unphysical divergencies as in [17]. The BCS self-consistency equation is recovered as an 
approximation from the flow. Neumayr and Metzner |23j study the interplay of d-wave 
superconductivity and Fermi surface deformations by a counterterm technique. The 
effective interactions causing the symmetry breaking are constructed order by order in 
perturbation theory. The counterterm is determined by setting the self-energy on the 
Fermi surface equal to zero, giving rise to a self-consistency equation. 

It is shown in [TT| [2^ that such self-consistency equations are also naturally 
contained in the fRG method applied here and in [TTl [T2| [23] . This method is exact 
for the models considered, as is the self-consistency treatment. However, it should be 
noted that for general models, self-consistent calculations just as fRG calculations are 
not exact anymore and their results may differ. It may depend on the specific case 
whether self-consistent calculations or the broader and more fundamental fRG method 
give the more useful results. 

The role of counterterms in the renormalization group literature has focussed on the 
removal of divergencies. In this work, the counterterm serves the purpose of selecting 
a certain symmetry-broken configuration. It also opens a gap in the spectrum of the 
system, which removes the divergencies on the Fermi surface treated in the literature. 

This paper is structured as follows. In section [21 we introduce the model, presenting 
and explaining the mean-field solution in section [3l We outline the renormalization- 
group scheme and introduce our conterterm extension in section |H where we furthermore 
provide a numerical RG analysis of the model from section [21 We discuss our work and 
provide an outlook in section [51 
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2. Model 

To illustrate the utility of the fermionic functional renormalization group (fRG) scheme 
introduced in section |U we consider a tight-binding model on a (i-dimensional hypercubic 
lattice with an interaction restricted to momentum transfers oi Q = (7r,7r,...). The 
Hamiltonian reads 



Hred = Cik) cicj. ct C-. -.ct C 

A^^^ ' k k AT Z-^ ki ki+Q k2 , 



]\f Z-^ ki ki+Q k2 k2-Q 
k k\,k2 

+ E (Ae.. + - A.) 4c,^^.. (1) 
k 

In this work, ^{k) = e{k) — fi includes a chemical potential /i, which can be adjusted to 
produce a system exhibiting a first-order phase transition. We assume in the following 
that the dispersion fulfills the nesting relation e{k) = —e{k — Q). The kinetic energy 
^{k), interaction, and external- field parts related to Agxt are also introduced and briefly 
discussed in Chapter 2 of [12]. For Vq > the interaction term can lead to charge- 
density- wave ordering with wavevector Q. This will be discussed below employing a 
mean-field approach. In fact, due to the restricted form of the interaction of the model, 
the mean-field treatment and the RG approach described below become exact in the 
thermodynamic limit, where the number of lattice sites N ^ oo (for a diagrammatic 
argument, see [21]). The same two-sublattice charge modulation is induced by the 
last term, which breaks the translational symmetry of the Hamiltonian by coupling an 
alternating charge field to the fermions. In the fRG treatment. Si will be used as the 
initial condition for the off-diagonal self-energy. The counterterm Ac will be included in 
the bare propagator to prevent this from affecting the physics of the system. We always 
set Si = Ac which guarantees a cancellation at the end of the fRG flow. Aoxt allows us 
to study the effect of an external field and hysteresis phenomena. 

3. Mean- field treatment 

The low-temperature state of the model without explicit symmetry breaking, i.e. for 
(Acxt + Si — Ac) = 0, depends on the chemical potential /i. For /i = and 0.5 particles 
(per spin orientation) per lattice site, the Fermi surface is perfectly nested and mean-field 
theory finds a second-order transition toward the charge- density wave state at a critical 
temperature Tc. In the ordered state, A = ^ Sfc(c^c^_,_Q) becomes nonzero. Similarly, 
the RG without self-energy corrections would find a runaway flow at a nonzero critical 
scale Ac. The runaway flow and the second-order mean- field transition are removed by 
a sufficiently large fi of the order of Tc{fi = 0). For such /x, a local stability analysis 
of the symmetric state via an expansion around A = does not detect any instability. 
However, below a transition temperature Tt, a mean-field search reveals global minima 
of the thermodynamic potential as function of A which do not develop out of the local 
minimum at A = when the temperature is lowered. As the temperature is lowered. 
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the system jumps from A = to the nonzero value, undergoing a first-order phase 
transition. The resuhing phase diagram is shown in Fig. [TJ The thermodynamic 
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Figure 1. /i-T phase diagram of ([T]) for Vq = 2.0i. 

potentials as functions of A for /i = and for first-order cases with /i 7^ are shown in 
Fig. [31 

To understand how such a simple model can have a first-order transition, we 
consider the self-consistency equation 

^ = ]^i:^r^^{m)-m)), (2) 

where f{x) = (exp(x/T) + 1)^^ is the usual Fermi distribution and the eigenvalues are 

E^ = ±R-fi, (3) 

where R = ■>/ (^fc)^ + A^. ([2]) is obtained from the differential of the grand canonical 
potential Q by demanding = 0. To be specific, we base the following discussion on 
a two-dimensional square lattice with nearest neighbor hopping t (which is used as the 
unit of energy for the remainder of the paper) at T = 0. The self-consistency equation 
is plotted in Fig. [2l The RHS has a kink at A = because one of the branches of 
the eigenvalues changes its sign on the whole umklapp surface for A = cutting 
out this part of the integration area of the RHS. It can be shown that at the kink the 
slope diverges as (1/ a/|/u| — A) log(|/i| — A) for A y At T = 0, the RHS does 
not depend on /i for A > since /i enters the RHS only via the Fermi functions. 
Since (E^) is always positive (negative) for A > the thermodynamically stable 
solution is always half filled in this case. The dependence of the filling on A should 
not be neglected, as is sometimes done in the literature [25]. In fact, mean-field theory 
already predicts the breakdown of the antiferromagnetic gap away from half filling in 
the homogeneous solution, in agreement with the quantum Monte-Carlo calculations 
of [25]. For a detailed discussion, see [H]. Since ([2]) is obtained from = 0, all 
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Figure 2. The right hand side of the gap equation ([2]) has a kink at A = leading 
to further solutions. The solutions at A = 0,0. 39t are minima while the solution at 
A = 0.24t is a maximum of Q.. All data obtained for Vb = 2.0i. 

intersection points are extrema (or saddle-points) of VL. If all intersection points are 
extrema of fi(A), the outermost one must be a minimum to guarantee the existence of 
a stable solution. By construction, the integral of the bisector minus the RHS of ([2]) is 
the energy gain through the opening of a finite gap. Thus, comparing the two shaded 
areas of Fig. [2] determines which of the two minima given by the intersections has lower 
energy. 

Plots of the thermodynamic potential are shown in Fig.O For small > 0, a local 
minimum develops at A = 0, which becomes the global and finally the only minimum 
if /i is increased further. Numerical study of ([2]) and VL for various T and leads to the 
phase diagram of Fig. [1] 

4. Renormalization group 

In this section, we outline an extension of the one-particle irreducible (IPI) fermionic 
functional renormalization group (fRG) scheme [Ij with the goal of studying first-order 
phase transitions. We start from the flow equations as suggested by Katanin [9j , studied 
for the case of the breaking of a continuous symmetry in |TT] and for the case of the 
breaking of a discrete symmetry in |12j. We use the interaction- flow method [16], which 
permits the fruitful employment of counterterms. We show flows for first-order and 
second-order phase transitions, discussing the effect of the counterterm on possible 
approximations. We also study the effect of an external field on the flow, identifying a 
criterion for unphysical starting points. 
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Figure 3. The thermodynamic potential for t' = 0, Vb = 2M. The position of the 
minimum at non-zero gap values does not change with /i, while the energy gain strongly 
depends on it. il(A = 0) is subtracted from all curves. 

4.1. Setup 

The basic setup of a IPI fRG calculation with Katanin's modification is outlined for 
two different models in [HI [T2]. To understand the problem solved in the following, 
we consider the order-parameter dependence of the thermodynamic potential in the 
plot for /i = in Fig. [3l In [H], [12], a small external field is included as the initial 
condition for the order parameter. This external field is relevant in Wilson's fixed- 
point picture, but does not appreciably change the physics away from critical scales. It 
biases the fRG calculation toward the thermodynamic potential minimum at finite A 
(see Fig. [3] /i = 0), which becomes the endpoint of the flow. This method is unable 
to deal with first-order phase transitions where the thermodynamically stable state is 
separated from the symmetric state by an energy barrier. A typical thermodynamic 
potential as function of the order parameter is depicted in the plot for /i = — 0.2t in 
Fig. [3l Because the external field must be small, the flow always veers toward the 
symmetric thermodynamic potential minimum. Other minima are inaccessible, even if 
their thermodynamic potentials are smaller. Furthermore, the method cannot be used 
to study hysteresis effects appearing if systems exhibiting second-order phase transitions 
(see the graphs for /i = and /i = — O.lt in Fig. [3l) are placed in an external field. The 
external field destroys the axis symmetry of f2(A), but always biases the flow toward 
the global minimum of the thermodynamic potential, making it impossible to study 
the metastable configuration. Thus, the challenge encountered here is to set up an 
fRG scheme with a parameter that allows the selection of any stable or metastable 
configuration as the endpoint of the flow. 

The basic idea that will be exploited in the following is to include in the calculation 
a counterterm Ac of arbitrary strength that cancels an equally strong external field. 
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but only at the end of the flow. In the fRG calculation, the external field is taken into 
account as the initial value Sj of the order parameter S, while the counterterm A,, is 
added to the naked propagator Qq. Schematically, the matrix Green's function of the 
system reads 

G-' = - E. (4) 

X 

where x is the cutoff function and Qo is the inverse bare propagator without counterterm. 
For any % ^ 1, counterterm and self-energy do not cancel each other. However, the 
special case x = 1 is approached at the end of the RG flow. Thus, the vertex function 
of the symmetric model is obtained. 

We first consider the case of x = XA(fc) describing the sharp momentum-shell cutoff 
of conventional RG schemes (see [1]). Then, only momenta for which Xa(^) = 1 are 
taken into account in the RG flow. For these momenta and E cancel exactly at the 
start of the flow, causing E to remain stationary. This is due to the right-hand side of 
the IPI flow equation for the order parameter being proportional to the anomalous part 
of the Green's function which vanishes if Ac and E cancel. This argument generalizes 
to the case of multiple order parameters, as can be seen by considering the adjugate 
matrix of the Green's function. Consequently, Ac and E also cancel at all stages of the 
flow. In the case of a second-order transition, the fRG will diverge at a nonzero scale 
since the spectrum remains unchanged. In the case of a first-order transition, the flow 
remains finite, signaling incorrectly the stability of the symmetric phase. We conclude 
that we have to employ a softer cutoff function. 

Next, we consider the softest possible cutoff function, provided by the interaction 
flow scheme [IE]. This scheme works by linearly turning on the interaction with the 
flow parameter (note the similarity to the treatment by Mahan [26]) while turning on 
the counterterm with the square root of the flow parameter. By rescaling the fields, 
this is found to be equivalent to employing a cutoff function x which is constant at all 
momenta, zero at the start of the flow, and one at the end of the flow. With this, the 
counterterm and the external field do not cancel in (jl]) until the end of the flow, and all 
configurations become accessible, as will be illustrated in the following. 

For a general cutoff function Xh with flow parameter A and a general model with 
two-particle bare interaction, the flow equations were obtained in [9]: 

S».-iT.(v^G»A(^)c^^) (5) 

= -- Tr fv^^ (G^G^) V^) + ©(V^), (6) 

where G^ = (Q/xa — E^) ^ and the traces run over all degrees of freedom of G^. The 
contribution of 0(V'^) to (jS]) vanishes in the case of the Hamiltonian ([T]) when taking 
the thermodynamic limit due to the restricted momentum structure of the interaction. 
For the interaction flow, the cutoff function is 



XA := VA, (7) 
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where A runs from zero to one. Considering e = ^ + yU as mapping the Brillouin zone 
onto M, denoting the effective gap Af := — Ac + x^? R '■= y^^^ + A^ , and := —fi±R, 
we obtain the A-dependent inverse free propagator 



[UJr-. 



1_ ( iun - e + n Ac 
X \ icun + e + /i 



(8) 



and the A-dependent Matsubara Green's function 

(9) 

We have omitted A-dependencies in the notation for brevity, as we continue to do in the 
following. ([H]) and allow the calculation of the single-scale propagator § = — GQG. 
Tracing over all degrees of freedom, analytically evaluating the Matsubara sum, and 
multiplying with the effective interaction and diagrammatic factors, we obtain the flow 
equation 




J (27r)''4v^l R [ 



Similarly treating the right-most diagram of Fig. 3 from [121 alternatively evaluating 
for the special case ([T]) yields the interaction flow equation 

^ ~ ^ J {27:Y2R^\ R ^ y^R ^ 

{f'{E~) + f'{E+)) (^3eA^ + A 

- (r (i^;-) + r (E+)) AA,2i?| . (11) 

We furthermore need to calculate the thermodynamic potential using the functional 
renormalization group. From equation (48) of [4J, using fii.a. = T70 and equations (39) 
and (45) of we determine the flow equation 

fi = |Tr((G-Q-i)( 

= Tr (1 + Q-^S + (Q-^S)' + ...)- ^ 



2 VX ^ ^ X 

-^Tr(s^G). (12) 



2 V X 

If the flow is started at the thermodynamic potential of the non-interacting system, the 
full thermodynamic potential is recovered. Since we will always do so, we have written 
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Vt instead of fii.a. in (1121) . The right-hand side of fll2p can again be drawn as a IPI 
diagram. Evaluating the Matsubara sums, f|T2l) reads 



( TTUj) . (fTTl) . and (fT^ constitute a closed system of integro-differential equations which 
can be numerically solved given appropriate initial conditions. 

4.2. Flows for first- order phase transitions 

We first consider a system exhibiting a first-order phase transition at a transition 
temperature Tt. Studying flows below Tt (see Fig. H]), we notice two strong attractors. 
By its lower thermodynamic potential at the end of the flow, one of them can be 
identified with the stable, symmetry-broken configuration. The values for the order 
parameter, effective interaction, and thermodynamic potential difference reproduce the 
exact mean-field results. Note that the final values do not depend on the magnitude 
of the counterterm. The dependence of the results on the external field, known from 
m EZl [m, is eliminated. 

Considering Fig. 111(a), we see that there is a separatrix between the effective gap 
flows to the metastable configuration and the flows to the stable configuration. For more 
complicated problems than ([1]) , the flow probably has to be stopped at A < 1 because of 
the rise of large effective interactions which are not correctly taken into account by our 
method. Therefore, flows which exhibit weak changes of the effective interaction away 
from A = 1 and weak changes of the order parameter close to A = 1 are most promising 
for determining order parameters in more complicated models. We see that the flows 
closest to the separatrix are not optimal flows since A (A) still shows a large slope at 
A = 1. This final slope is smallest for a counterterm of roughly twice the magnitude of 
the order parameter. Such a flow would yield an excellent approximation of the order 
parameter even if stopped at A = 0.75. 

Considering the flows of the effective interaction in diagram |4](b), we do not find 
a separatrix as in (a). Instead, the flows corresponding to the ones close to the 
separatrix in (a) develop a shoulder which becomes a maximum for certain values 
of the counterterm (for an example of such a maximum, see diagram [5](b)). The 
flows exhibiting the strongest effective interactions correspond to the flows giving the 
worst approximations for the order parameter if terminated prematurely. This entails 
for applications to models where this method is not exact that carefully chosing the 
counterterm can significantly simplify the calculation and improve the accuracy of the 
results. 

The flows of the thermodynamic potential shown in Wic) again exhibit separatrix 
behaviour. This is in contrast to the flows above Tt. For the thermodynamic 
potential, the flows closest to the separatrix yield the best approximations if terminated 
prematurely. However, it is apparent from|l](c) that the flow must be continued until the 
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Figure 4. Flows for Vq = 2t at ^ = 0.245t, T = O.Olt < Tt, and Ac increasing from 
O.Olt to l.Olt in increments of O.li. Broken lines denote flows converging to the stable 
symmetry-broken configuration; solid lines denote flows converging to the metastable 
symmetric configuration. The flow starts at A = and finishes at A = 1. (a) Effective 
gap Af. Ac for each graph can be read off at the y-axis. (b) Effective interaction V . 
Flows converging at 10.14t pass through greater values with increasing counter-term. 
Flows converging at 4.43t behave inversely, (c) Thermodynamic potential Vl. Flows 
pass through smaller values with increasing counterterm. 



external field is completely compensated by the counterterm to obtain reliable results 
for the thermodynamic potential. 

Flows above Tt are illustrated in Fig. [51 Again, we clearly discern two attractors. In 
this case, however, the thermodynamic potential flows of the symmetry-broken-phase 
attractor cross the thermodynamic potential flows of the symmetric-phase attractor. 
The symmetric phase is therefore thermodynamically more stable. Apart from this, the 
flows behave similarly as those below Tt. For inconveniently chosen counterterms, the 
maxima in the effective- interaction flows are clearly visible in|5]^b). Nevertheless, the 
flows reproduce the exact mean-field results. 
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Figure 5. Flows for Vq = 2t at /i = 0.245t, T = 0.04t > Tt, and increasing 
from O.Oli to l.Oli in increments of O.lt. Broken lines denote flows converging to 
the metastable symmetry-broken configuration; solid lines denote flows converging to 
the stable symmetric configuration. The fiow starts at A = and finishes at A = 1. 
(a) Effective gap Af. Ac for each graph can be read off at the y-axis. (b) Effective 
interaction V . Flows converging at 11.90t pass through greater values with increasing 
counter-term. Flows converging at 5.62t behave inversely. (c) Thermodynamic 
potential O. Symmetric flows pass through smaller values with increasing counterterm, 
and so do symmetry-broken flows. 



J^..3. Flows for second-order phase transitions 

The flows for second-order pliase transitions beliave similarly to the flows for first-order 
phase transitions, but only a single attractor appears, as can be seen in Fig. [6](a). Fig. 
Mjo) shows that it is possible to suppress the effective interaction during the flow by 
chosing a large counterterm. Approximations for the effective gap can be obtained by 
stopping the flow before A reaches 1. The quality of such an approximation depends on 
the counterterm chosen. Again, the best approximation can be obtained by chosing a 
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Figure 6. Flows for Vb 2t at fi=^Q,T « Tc, and Ac increasing from O.Olt to l.Olt 
in increments of O.lt. Broken lines denote flows converging to the stable symmetry- 
broken configuration. The flow starts at A = and finishes at A = 1. (a) Effective 
gap Af. Ac for each graph can be read off at the y-axis. (b) Effective interaction V. 
Flows pass through smaller values with increasing counterterm. 



counterterm of twice the non- approximated value of the gap. 
4-. 4- External field 

If the external field Aext in dl]) is zero, the half-filled system below exhibits 
two degenerate stable configurations (local minima of the thermodynamic potential) 
distinguished by the sign of the order parameter. If Aext is non-zero, this degeneracy 
is lifted. The fRG scheme outlined above allows us to select the endpoint of the flow 
independently of the external field by appropriately setting the counterterm, in contrast 
to the scheme from p!T| [12] . The flows toward the metastable configuration are shown 
in Fig. [7] (broken lines). As the external field is increased toward a critical value, the 
order parameter for the metastable configuration vanishes and the corresponding local 
minimum of the thermodynamic potential disappears. Beyond this critical value, the 
flows for negative counterterms become divergent as the effective gap reaches values 
close to zero, exposing low-energy modes. Such a divergence can be used as an indicator 
for an adversely chosen starting point. The flows for positive counterterms (full lines in 
Fig. [7]) always attain the stable solution. Calculating the flows as for Fig. [71 but for 
negative values of the external field Aext, we obtain the hysteresis curve of Fig. [HI 

5. Discussion and outlook 

We have described a counter-term scheme within the framework of the functional 
renormalization group concept for interacting fermions. The scheme is able to detect 
ordered ground states which, in terms of the effective bosonic potential for the ordering 
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Figure 7. Flows of Af at ^ = for Vq = 2t. Full (broken) lines have Ac = (-)i. The 
external field increases from to 0.15i with increasing final values of the flow. The 
flow starts at A = and finishes at A = 1. 
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Figure 8. Hysteresis effect in an external field. All solutions represent symmetry- 
broken configurations. Green lines (x -symbols) represent stable, red lines (-1— symbols) 
mctastable configurations. 



field, occur as first-order transitions. Second-order transitions are captured as well. 
Therefore, the new scheme remedies one shortcoming of conventional functional RG 
techniques (e.g., [28l [29t [30] ) for interacting fermions which are blind with respect to 
symmetry-broken states which are separated by an energy barrier from the symmetric 
state. A second shortcoming that is addressed is the need for a finite external field 
p!T| [T2] . which does not arise in our approach. As a third improvement, at least 
for discrete-symmetry breaking away from the critical temperature, our scheme can 
circumvent the development of large interactions in the flow. 

The validity of the new scheme has been demonstrated in a mean-field model, 
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where our truncated version of the fRG equations and mean-field theory are known to 
be exact. More generally, the scheme can be used to target individual minima of the 
effective bosonic potential, for example to study metastable states in the presence of an 
external field. As an example, we have shown the calculation of a hysteresis curve. 

The fRG scheme can be applied to more realistic models such as the Hubbard model 
as well. Here a new complication arises, as the truncated form of the RG equations 
used is no longer exact. Hence large values of the coupling constants could spoil the 
validity. This is a general problem for perturbative functional renormalization group 
schemes which can for example be circumvented by partial bosonization at the expense of 
introducing a dynamical bosonic field [El [71 El El] . The counter-term technique presented 
here allows to avoid this problem in certain cases. As demonstrated and described above, 
in the case of the breaking of a discrete symmetry, no large couplings occur unless one is 
close to a curvature change of the effective bosonic potential. Hence for these cases, an 
inclusion of non-mean-field-type interactions should be possible without much trouble. 
For broken continuous symmetries, Goldstone modes may again lead to large couplings. 
It will be interesting to see whether the Goldstone physics only enters the flow in a 
critical way toward the end of the flow, when the off-diagonal selfenergy has already 
converged to the final value. The results above indicate a rather quick convergence for 
most values of the counter-term. These issues will be the subject of future work. 

We emphasize however that even in the case that the flow cannot be followed over 
the full range, the present scheme should be useful for the detection of symmetry-broken 
states which are separated by an energy barrier from the symmetric state. For example, 
if the conventional analysis signals a locally stable symmetric state, the counter-term 
technique can be used to scan non-zero values of the suspected order parameter. If 
there is a lower thermodynamic potential minimum further out, the self-energy shows 
a flow toward non-zero values. Further, there should be a local curvature change of 
the effective bosonic potential somewhere between and the value related to the global 
minimum, which is detected as a runaway flow of the interactions in the corresponding 
channel. Hence, the method is useful at least for qualitative scans of the bosonic 
potential landscape. More generally, variations of the proposed scheme can be used 
for the perturbative analysis around states which differ from the free ground state of 
the model. This as well should lead to numerous new applications. 
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